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Abstract. - The interplay between self-difi^usion and excitation lines in space-time was re- 
cently studied in kinetically constrained models to explain the breakdown of the Stokes-Einstein 
law in supercooled liquids. Here, we further examine this interplay and its manifestation in 
incoherent scattering functions. In particular, we establish a dynamic length scale below which 
Fickian diffusion breaks down, as is observed in experiments and simulations. We describe 
the temperature dependence of this length scale in liquids of various fragilities, and provide 
analytical estimates for the van Hove and self-intermediate scattering functions. 

A ten-day journey starts with a single step. — Laotse, Tao Te King. 

In this paper, we consider the process of self-diffusion of probe molecules in supercooled 
liquids. Figurenshows the trajectory of a such a probe obtained from a model of a supercooled 
liquid [1]. At conditions shown, the structural relaxation time of the model is of the order 
of 10^ microscopic time steps. The left panel of Fig. ^ extends over this range of time. The 
right panel extends three orders of magnitude longer in time, and one order of magnitude 
larger in space. Here, the trajectory looks like a random walk of Fickian diffusion, unlike 
the trajectory in the left panel. In this paper, we describe the crossover from non-Fickian to 
Fickian diffusion, and identify the length scale, , that characterizes the crossover. 

In our perspective supercooled liquids are modeled by a coarse-grained mobility field evolv- 
ing with simple empirical rules [2] . An essential prediction of our modeling is that dynamics 
becomes spatially correlated [3], i.e., the growth of timescales is accompanied by the growth 
of dynamical lengthscales, giving rise to the phenomenon of dynamic heterogeneity observed 
in experiments and simulations [4-7]. This finding suggests a degree of universality and thus 
utility for a coarse-grained perspective. Nevertheless, it is important to test of our approach by 
using it to revisit all sorts of experimental and numerical studies of supercooled liquids [1,8-10]. 
Generic properties are only weakly dependent upon details of the models, which become im- 
portant, however, for quantitative comparisons to experiments or simulations [8]. In most of 
this paper we will therefore pursue our investigations in the simplest lattice model of this fam- 
ily [2], namely the one-dimensional Fredrickson- Andersen model (hereafter ID FA model) [11]. 
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Fig. 1 - A typical trajectory of a probe molecule in the ID FA model at T = 0.25 is shown twice on 
different time and length scales. The a-relaxation time is ~ 10^, the diffusion constant D ~ 10~*. 



This is defined by the Hamiltonian H = rii, where rii = 0, 1, with dynamics constrained by 
isotropic dynamic facihtation. A site needs at least one neighbour with = 1 to change state, 
with Boltzmann probability. The concentration of excited sites is c = {ui) = [1 -I- e^^"'"]^^, T 
being the temperature. Lengths are in units of the lattice spacing which we set to one. 

Self-diffusion is studied by introducing a probe molecule [1], which makes it possible to 
map 'spins' back to 'particles'. The probe's position, x{t), evolves with reduced time as 
x{t -|- 1) = x{t) ± nx{t)nx±i{t), i.e., the probe is allowed to move only when sitting on an 
excitation and jumps to a neighbouring excited site. Due to coarse-graining, fast vibrations are 
ignored. In that sense, our results are close to dynamic studies of inherent structures [13,14], 
but this is not crucial as far as the long-time relaxation is concerned. Numerically, we studied 
the dynamics of probe molecules in the ID FA model using Monte Carlo simulations. At each 
temperature, the system size, L, is chosen to be much larger than the mean distance between 
excitations, £(T) = l/c(T), in order to avoid finite size effects. We typically used L > 10£(T) 
and averaged over at least 10^ independent probe trajectories. We have also studied the 
behavior of a probe in the East model [12], a more constrained and fragile counterpart of the 
ID FA model, and results are briefly discussed near the end of the paper. 

The trajectories in Fig.^show the typical time evolution of the position of a probe molecule 
in the ID FA model at a fixed, low temperature. Due to the kinetic constraint, excitations in 
dynamically facilitated models propagate as continuous excitation lines in space and time [1,3]. 
On the a-relaxation timescale motion is non-Fickian with long periods where the molecule is 
immobile, seen as plateaux in Fig.^ punctuated by shorter periods where it travels over a few 
sites. Plateaux correspond to the particle being 'caged' in regions of space far from excitations, 
and having to wait for an excitation line to reach its position and allow it to move. As discussed 
in Ref. [1], the average time it takes a probe to be reached by an excitation line and thus move 
for the first time is the mean persistence time, which is also the a-relaxation time, r^, of 
these models. Further, once the particle has moved, subsequent steps occur when it is hit 
again by excitation lines, i.e., each time there is a microscopic exchange event. The average 
exchange time, t^, determines the self-diffusion constant, D. Since grows much more 
slowly with decreasing temperature than Tq, the breakdown of the Stokes-Einstein relation 
follows [1]. This means also that the propagation of the particle is given by a competition 
between these two fundamental processes: persistence and exchange. Which one dominates 
at a given temperature will depend on the time and length scales of motion. 

To quantify the above observations, we study the histogram of particle displacements, i.e. 
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Fig. 2 - Left: Ga{r,t), Eq. at fixed temperature and various times has a bimodal form from 
slow/fast molecules. Right: Diffusive part of Gs{r, t) rescaled after Eq. ^ shows decreasing deviations 
from the Gaussian (dashed line) with increasing time. 



the selLpart of the van Hove function, 

Gs{r,t)^(6{r~[x{t)-xm)), (1) 

where brackets indicate an average over trajectories. The above discussion suggests that at 
times of the order of Tq,, some particles are still trapped in their initial position, while oth- 
ers have already escaped and rapidly diffused over a certain distance. The distribution (Q) 
should reflect this bimodal character with a contribution in r = from the 'slow', immobile 
particles, and contributions at non-zero displacements from diffusing, 'fast' particles. Numer- 
ical results in Fig. |21 (left panel) fully confirm this expectation. Coexistence of fast and slow 
subpopulations of particles is a major indicator of dynamic heterogeneity [4-7, 15-17]. 

The behavior of the van Hove function can be anticipated from a simple analysis. Each 
time the probe particle is allowed to move it makes a random walk step. Therefore, it performs 
a random walk, but the time lag between steps fluctuates. The time lags are determined by 
the dynamics of the mobility field of the host liquid. For ease of analysis, we neglect the back 
reaction of the probe on this field. In this approximation, the probe process is a random walk 
at random times, or continuous time random walk [18]. The probability for the particle to be 
at position r at time t is then Gs{r,t) = Y^^=o''^rn{t)(l)^"^\r), where iTm{t) is the probability 
that the probe jumped m times up to time t, and (j)^™^\r) is the probability that a random 
walker is at site r after m steps, having started from the origin. 

The first step is a persistence event. If p{t) is the distribution of persistence times [9], 
then the probability 7ro(t) of not having an event up to time t is given by 7ro(t) — P{t), where 
P{t) = dt'p{t') is the persistence function [3,9]. If we further assume that successive 
exchange events are uncorrelated, which is a good approximation in the FA model [1], then 
■Km{t) can be written as a multiple convolution of the exchange time distribution ^(f) [19]. The 
Laplace transform of -Kmii) for m > then reads: TTm{s) = p(s)i/''"~-^(s)[l — ■0(s)]/s, where 
hats indicate Laplace transforms. If we also Fourier transform in space we obtain Fs{k,s), 
which is the Laplace transform of the self-intermediate scattering function, Fs{k,t), and reads 

Fs{k,s)=P{s)+ cos {k)^ L^ilf^. (2) 

s 1 — cos (fc)?/;(s) 



This is the MontroU- Weiss equation for the propagator of the probe particle [18]. From Eq. ^ 
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we immediately obtain a k and T dependent timescale: T{k,T) — lims^o Fs{k, s). Since 
the distributions p{t) and iplt) are both narrow (all moments are finite [3]), their Laplace 
transforms read, for small s, p(s) rs 1 — TaS and 'ip{s) « 1 — r^s, which leads to: 

For small k this gives r w + (k^D)^^, since D = 1/{2tx) [1]. 

The inverse Laplace transform of the factor cos (fc)[l — V'] [s — s cos (fc)'?/']"^ in Eq. Q defines 
a propagator for diffusing molecules, i<diff(fc,t). Due to the narrowness of tp, it becomes that 
of normal diffusion in the large t, small k limit, _Fdiff(fc,i) ~ ex.p{—Dk^t). Inverse Laplace 
transforming ||2Jl we obtain Fs{k, t) — P(t) + dt'p(t')Fditf{k, t — t'). In the regime of strong 
decoupling this can be further approximated by: 

Fs{k,t) ^ P{t) + [\ ~ P{t)]FMk,t) . (4) 

The corresponding van Hove function then reads, 

G.(r,i) « P{t)5{r) + i-£^ exp f-^) , (5) 



where we have approximated G'diff(?', i) by a Gaussian for all times. We find that Eq. (jSJ 
works well for all times and temperatures, see Fig. |21 (right). Deviations from the Gaussian 
tend to disappear at long times, as expected. The mean square displacement reads (r^) = 
dr r^Gs{r,t) ~ [l — P{t)^a{t), where a{t) is the mean square displacement obtained by 
restricting the average to the fast particles. For large times, a{t cx)) « (r^) « 2Dt, while 
a{t ^ 0) = 1 by construction. Physically, Eq. (01 shows that it is a good approximation to 
think of molecules as existing in one of two subpopulations of immobile and mobile particles. 

Our primary results in this paper all follow from Eq. It is valid only for lengths and 
times larger than our lattice spacing and unit of time, respectively. Thus the delta-function 
in Eq. (jSJ contributes to an experimentally observed Gs (r, t) only after coarse-graining and 
would be replaced at smaller length scales by a smoother distribution, typically Gaussian. 
A two-Gaussian fit to the van Hove function is reported in confocal microscopy experiments 
performed with colloidal suspensions [15]. If one is predisposed to think of the distribution in 
terms of a single Gaussian, it will appear as if the distribution has 'fat tails.' Such tails have 
been described as indicators of dynamic heterogeneity [5, 16, 17]. 

The approximation JSJ also provides a simple explanation for the time dependence of the 
non-Gaussian parameter, a2{t) = ■^{r'^) / {r'^) ■ By construction, a2{t) = 1 for a Gaussian 
Gs (r, t) . Instead, 02 (t) reflects the transition from a Gaussian (fast vibrations at short times) 
to another Gaussian (Fickian diffusion at long times), the distribution being non-Gaussian 
at intermediate times when both terms in ^ compete. From Eq. we obtain a2{t) = 
[1 — P{t)\ , which is monotonically decreasing with time. This expression is valid beyond a 
coarse graining time, and therefore does not describe the very short-time Gaussian behaviour 
that would make a2(t) non-monotonic. A similar monotonic behaviour is found in molecular 
dynamics studies of inherent structures dynamics [14]. This expression does, however, explain 
why a2 (t) can reach values much larger than unity. 

Figure 01 presents our numerical results at various temperatures and wavevectors for the 
self-intermediate scattering function Fs{k,t). At high temperature, T = 1.5, relaxation is 
homogeneous and Ff,{k,t) decays exponentially at all wavevectors. The effect of dynamic 
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Fig. 3 - Self-intermediate scattering function in the ID FA model, for different T and the same range 
of fc: k = 27r/n, with n = 2, 4, 6, 8, 10, 15, 20, 30, 40, 60, 80, and 100 (from left to right). Dashed 
lines at r = 0.3 represent the persistence P{t) and the diffusive behaviour exp [-{2Tr/80fDt]. 

heterogeneity becomes fully visible when T decreases. First, time decays are non-exponential 
at large wavevector. Second, Fs{k,t) becomes exponential when k decreases at constant T, 
as suggested by our observations of Fig. ^ Finally, curves for different wavevectors that are 
distinct at high T superpose when T is decreased, and Fs{k,t) becomes ^-independent. 

These observations are rationalized by Eq. Q). For large fc, the last term in Q is expo- 
nentially suppressed, and Fs{k,t) « P{t), which is exp(— -^/t/rc) in the ID FA model [20]. 
This behaviour justifies our identification of Tq with the persistence time [1,9]. At small k 
and large times, Eq. becomes Fs{k,t) « exp(— fc^Dt), as confirmed by Fig. O At in- 
termediate fc, mixed behaviour between diffusion and persistence is predicted, and observed. 
These observations naturally imply that stretching is also a fc and T dependent property, at 
odds with predictions stemming from mode-coupling theory [22] , which indeed poorly describe 
self-dynamics at intermediate wavevectors [22]. 

The crossover from persistent to Fickian diffusive dynamics shows up in the behaviour of 
the relaxation times T(fc, T), as given by Eq. Q. Numerically, we extract the timescales from 
the scattering functions in the usual way, Fs(k,T) = e^^. Results are presented in Fig. ^ 
Fickian behaviour, r w k~^D~^, is obeyed at all fc at high T. When T decreases, diffusive 
behaviour is restricted to smaller and smaller fc, and is progressively replaced, at larger fc, by 
the fc-independent form r « Tq,. This leads us to define a temperature dependent characteristic 
length scale, i*{T), which determines the onset of Fickian diffusion. From Eq. © we get 

t ~ /d^ . (6) 

It is therefore possible to collapse all timescales on a master curve rescaling times by k^^D^^ 
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Fig. 4 - Left: Relaxation times in the ID FA model, as a function of wavevector for different temper- 
atures, T = 1.5, 0.8, 0.6, 0.3, 0.2, and 0.15 (from bottom to top). The straight lines correspond to 
T ~ and T ^ const, the limiting behaviours predicted by Eq. I^. Right: timescales rescaled by 
the diffusive limit k~^D~^ are collapsed using the scaled variable k£* for both ID FA (filled symbols) 
and Ecist (open symbols) models. Data for the East model cover about 10 decades in timescales, with 
T € [0.31, 0.8] and the scaling curve has been shifted by a factor 5, for clarity. 



and space by £* . This is shown in the right panel of Fig.^ where £* is estimated from Eq. 
In Fig. ^ we also show a similar collapse of time scales that we have found when the East 
model [12] is used in place of the ID FA model. The East model is the fragile counterpart of the 
ID FA model. We see that the scaling behaviour predicted by Eqs. (OEl works well also in this 
case, despite the fact that in the East model successive exchange events are not uncorrelated 
[1], which was one of the assumptions in the derivation of Eqs. lO-lO- A behaviour similar 
to that of Figs. 13 and 21 was recently reported in molecular dynamics simulations of a binary 
Lennard- Jones mixture [23] , and in experiments on supercooled TNB [24] . 

For both the FA and East models the diffusion constant obeys a fractional Stokes-Einstein 
law, D ~ r-«, with ^ < 1 [1]. In the case of the FA model, ^fa = 2/A « 2/3, 2/2.3, 2/2.1 
for dimensions d = 1, 2, 3, respectively, A being the time exponent ^ c~'^ [21,25]. The 
Stokes-Einstein law, ^ = 1, is recovered for d > 4, the upper critical dimension of the FA 
model [21,25]. For the East model numerical results indicate ^East ~ 0.7 — 0.8, independent 
of d up to the highest dimensionality studied, d < 6. Equation (jSJ can be rewritten, 

t ^ r(i-«)/2 . (7) 

Therefore, £* will diverge when T ^ (or c ^ 0) if there is Stokes-Einstein breakdown, i.e., 
when f < 1. The Fickian crossover length, i*, measures something related to but distinct 
from the largest dynamic heterogeneity length, £{T), which can be measured for example 
through multi-point dynamic structure factors [6]. For both strong (FA) and fragile (East) 
systems, this length goes as £ ^ c^'^ with a spatial exponent v [3,25], and therefore always 
diverges as T — > in an Arrhenius manner. In general, we have that £ ^ £* , and in particular 
for fragile systems, Eq. shows that £* will grow faster than £, in a super- Arrhenius way. 
The reverse is true in strong systems where £ should grow faster than £* . In Ref. [23], the 
typical length scale £ of dynamic heterogeneity was used to rescale wavevectors in the analog 
of Fig. 0] Present results show instead that £ and £* are different quantities, although they 
might be hard to distinguish on a restricted temperature window. Finally, £* was identified in 
Ref. [26] by dimensional analysis which confuses i and £*. This particular confusion is avoided 
in Ref. [27] where the crossover length is also obtained, by assuming a memory function is 
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analytic for small wave vectors and its Taylor series can be truncated at second order even 
for wave- vectors that are not small. 

The physical basis for the Fickian crossover is that persistence time dominates over ex- 
change time. The former is the time for the first dynamical step, while the latter is the typical 
time scale of subsequent motion. This mechanism naturally leads to the breakdown of the 
Stokes-Einstein law [1], non-Gaussian van Hove functions with large tails, and sub-populations 
of fast and slow particles. The present picture is different from the idea that relaxation and 
diffusion result from averaging a time and its inverse, respectively, over a single distribution of 
relaxation times [27,28]. Our explanation that physically distinct processes (persistence and 
exchange) compete is consistent with the experimental observation [24] of a strong decoupling 
in a material with self-similar distributions of relaxation times. 
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